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ABSTRACT 

We present measurements of the two-point galaxy angular correlation function w(6) in the COSMOS 
field. Independent determinations of w(0) as a function of magnitude limit arc presented for both the 
HST ACS catalog and also for the ground-based data from Subaru and the CFHT. Despite having 
significantly different masks, these three determinations agree well. At bright magnitudes (Jab < 22), 
our data generally match very well with existing measurement s and w ith mock catalogs based on semi- 
analytic galaxy formation calculations of[Kitzbichler & White (2006) from the Millennium Simulation. 
The exception is that our result is at the upper end of the expected cosmic variance scatter for 6 > 10 
arcmin, which we attribute to a particularly rich structure known to exist at z ~ 0.8. For fainter 
samples, however, the level of clustering is somewhat higher than reported by some previous studies: 
in all three catalogues we find w(6 = V) ~ 0.014 at a median Tab magnitude of 24. At these very 
faintest magnitudes, our measurements agree well with the latest determinations from the Canada- 
France Legacy Survey. This level of clustering is approximately double what is predicted by the 
semi-analytic catalogs (at all angles). The semi-analytic results allow an estimate of cosmic variance, 
which is too small to account for the discrepancy. We therefore conclude that the mean amplitude of 
clustering at this level is higher than previously estimated. 

Subject headings: cosmology: observations — cosmology: large scale structure of universe — cosmol- 
ogy: dark matter — galaxies: formation — galaxies: evolution — surveys 
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1. INTRODUCTION 

The COSMOS field dScoville et al.ll27)07ch is the largest 
contiguous multi-wavelength probe of the high-redshift 
galaxy distribution, and a major task for the survey will 
be to extract improved measurements of galaxy cluster- 
ing at these early times. In this initial paper, we will be 
concerned with the simplest of these measures: the angu- 
lar two-point correlation function, w(0). Demonstrating 
a robust measurement of this quantity is a minimum re- 
quirement for verifying that the survey completeness is 
understood, as a basis for future more elaborate analysis. 
The main aim of this paper is therefore to present mea- 
surements of the two-point galaxy clustering statistic on 
the COSMOS field using three independently generated 
catalogs, and to compare the results with existing data. 
We also compare the amplitudes we measure to those 
found in a semi-analytic model of galaxy formation. 
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The key feature of the COSMOS field is that it is 
completely covered by the largest existing mosaic of im- 
age tiles from the Advanced Cam era for Surveys (ACS ) 
on the Hubble Space Telescope (jScoville et al.l l2007al) . 
With respect to ground-based surveys, in addition to the 
exceptional image quality, the great advantage of this 
dataset is the superior photometric accuracy and sta- 
bility over the entire field of view of the survey, which 
in turn makes it possible to measure the clustering of 
galaxies on large scales and at faint magnitudes where 
the amplitude of the galaxy correlation function w(9) is 
very small. The COSMOS field is currently the survey 
that probes the largest comoving scales at redshifts of 
around one. The total area covered, 2 deg 2 , does not ex- 
ceed existing studi es , in particular the UH8k study (1.5 
deg 2 . I Wilson! I2003T) and the shall ower DEEP2 measure- 
ments (5.0 dee 2 . iCoil et"aT][200l . However, COSMOS 
offers the unique combination of large contiguous area 
and depth, and also has the virtue of several indepen- 
dent and quite different imaging datasets in the same 
field. 

We use this rich dataset to investigate the clustering 
properties of the field galaxy population on degree scales. 
In future papers, we will present a more detailed study of 
galaxy clustering using photometric redshifts which can 
be used, for example, to divide our galaxy catalogues 
by type and apparent magnitude. Our objective here 
is simply to present the global properties of the field in 
terms of simple two-point statistics for catalogs selected 
by apparent magnitude. As our aim is to demonstrate 
the robustness of the results, we restrict ourselves to the 
i-band data, which is available for all three datasets con- 
sidered here. 

The COSMOS field has been imaged by many ground- 
based facilities , amon gst them the Subaru telescope 
(jTaniguchi et al.ll2007h . an d we compare the COS MOS 
ACS catalog (described in iLeauthaud et al.l I2007T ) with 
two ground- based catalog s : the Subaru optical catalog 
described bv lCapak et all (|200l and the CFHTLS-T03 
catalog used by McCracken et al. (2007, in prepara- 
tion). Each of these have significantly different masks, 
especially with regard to ghosting around bright stars. 
As will be demonstrated, the results of these indepen- 
dent determinations are in good agreement, and display 
a consistently higher amplitude at faint magnitudes than 
has been suggested in previous work. 

2. CATALOGS AND METHODS 

All of the three catalogs were prepar ed independently. 
For fu ll details of the ACS catalog, see ILeauthaud et al. 
( 20071) : the Subaru catalogs are described in lCapak et al. 
(|2007f i. A full description of the CFHTLS catalogs, based 
on the images corresponding to the release CFHTLS- 
T03, can be found in McCracken et al 2007. 

The ACS data are constructed from a mosaic of 575 
image tiles taken over 588 orbits of the HST. The 50% 
completeness limit of the catalog is 26.6 F814W magni- 
tudes. 

SuprimeCam consists of 10 Lincoln-Labs 8k x 4k 
CCDs with a plate scale of 0.2" per pixel. However, 
a plate scale of 0.15" per pixel was used for the final 
image to ensure that the images with good seeing were 
not undersampled. A special dither pattern including 
camera rotations (Subaru is an alt-azimuth telescope) 



was used to ensure every portion of the field was imaged 
by at le ast four different CCD s. Further details can be 
found in lTaniguchi et al.l(|2007l ). The Subaru catalog was 
based on a mosaic of 115 Subaru images taken with the 
Suprimc camera on the Subaru 8m telescope. The image 
has a median seeing of 0.6" and a 50% completeness of 
i = 27 A. 

In order to produce the best possible catalogue for cor- 
relation function measurements, we re-extracted a cata- 
logue ourselves from the Subaru tiles produced by Capak 
et al. We downloaded each tile and assembled them into 
a single large mosaic using the TERAPIX software tool 
swarp; we carried out the same procedure for the RMS 
maps. Following this, we used sextractor to extract a 
catalogue. Star-galaxy separation to i = 21.5 was per- 
formed using the f luxjradius compactness parameter 
which measures the radius which encloses 50% of an ob- 
ject's flux. Bright stars and defects were masked on the 
images. 

The CFHTLS catalog is derived from the TERAPIX 
CFHTLS-T03 release. The CFHTLS stacks were taken 
using the Megacam camera on the 3.6m CFHT telescope. 
Megacam covers 1 deg 2 with 0.205" pixels using 36 sep- 
arate 2048 x 4096 Rockwell CCDs. Note that, unlike the 
Subaru and ACS data, the CFHTLS-T03 image consists 
of a single Megacam pointing. The CFHTLS-T03 D2-i 
dataset comprises 153 images and has a median seeing of 
0.9". The 50% completeness of this dataset is Iab — 25.7 
magnitudes 22 . 

In all three catalogs, the star-galaxy separation is car- 
ried out using a morphological classifier. We empha- 
size that the classifiers were determined independently 
for each catalog; we did not, for example, use the ACS 
morphologies to perform star-galaxy separation on the 
CFHTLS or Subaru images, or use the Subaru images as 
detection images for the ACS data. The Subaru images 
cover the full 2 deg 2 of the COSMOS field while the ACS 
tiling covers a total of 1.7 deg 2 ; the CFHTLS-T03 images 
cover just the central 1 deg 2 of the field. 

In each catalog, regions around bright stars and near 
the edges of the field were masked. For the ACS cata- 
logs, we use the same set o f masks that were used for 
weak lensing measurements (jMassev et al.ll200"7T ). These 
masks also remove many blended objects. After mask- 
ing, in the Subaru catalogue there are 134,397 galaxies 
in the magnitude range 20 < i < 25; in the ACS cat- 
alogue and there are there are 124,665 galaxies in the 
same interval in magnitude. For the CFHTLS, there are 
52,521 galaxies in the magnitude range 20 < i < 24. The 
effective areas (total available area after masking) of the 
three surveys (Subaru, ACS and CFHT) is 1.6, 1.5, and 
0.7 deg 2 respectively (i.e. completeness of 80%, 88% and 
70%). We have experimented with varying the degree of 
masking by 'growing' the mask to eliminate pixels adja- 
cent to masked pixels. The results are robust even when 
> 50% of the area is masked. 

Figure [T] shows the galaxy number counts extracted 
from the three catalogues. The dotted lines indicate 
the magnitude limits adopted in this paper. The slight 
'break' in the counts at i* ~ 21 is an artefact caused by 
our morphologically based star-galaxy separation. The 
ACS count are slightly lower at i ~ 22 as a consequence 

22 For details see http: //terapix. iap. f r/cplt/tab_t03ym.html 
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of the improved star-galaxy separation in this dataset. 



3. CLUSTERING MEASUREMENTS 

We selected galaxies in progressively fainter slices of 
apparent i magnitude. For the purposes of the paper 
we assume that the instrumental AB total magnitudes 
measured in each catalog are eq uivalent; this is approxi- 
mately true. ICapak et al.l (|2007f ) present a detailed com- 
parison of galaxy photometry between the three catalogs 
described here. Figure 8 of their paper demonstrates that 
total instrumental magnitudes in each catalog agree well, 
to within 0.05 magnitudes. 

For each slice, we measure w for at range of angular 
separations 9 to 9 + 59 in a serie s of logarithmica ll y sep - 
arated bins using the standard lLandv fc Szalavl (|1993l ) 
estimator, 

... DD-2DR + RR 

W{9) = RR (1) 

with the DD, DR and RR terms referring to the number 
of data-data, data-random and random-random galaxy 
pairs between 9 and 9 + 59. The fitted amplitudes 
quoted in this paper assume a power-law slope for the 
galaxy correlation function, w{9) = A w (9 /deg)~ 5 ; how- 
ever this amplitude must be adjusted for the 'integral 
constraint' correction, arising from the need to estimate 
the mean galaxy den sity from the sample itse lf. This can 
be estimated as (e.g. [Adclbcrgcr ct al. 20051). 



C =^2 1 1 "'(0>'/<>i 



(2) 



where f2 is the area subtended by each of our survey 
fields. For the COSMOS field, We find C ~ 1A W by 
numerical integration of Equation [2] over our field geom- 
etry and assuming that galaxies closer than 1" cannot be 
distinguished. 

We used a sorted linked list in order to reduce com- 
puting time given the very large number of objects in 
each slice. These results are compared in Figure [3] The 
solid line shows measurements from the ACS; the trian- 
gles and stars correspond to measurements from Subaru 
and CFHT. At each angular bin for each survey, the er- 
ror bars plotted are simple bootstrap errors. Although 
these are not in general a perfect substitute for a full 
estimate of cosmic variance (e.g. using an ensemble of 
simulations), th ey should give the correct magnitude of 
the uncertainty (|Mo et al.lll992l ). In particular, these au- 
thors show that bootstrap errors yield sensible uncertain- 
ties on power-law fits to correlation-function data when 
the points are treated as independent. 

For each of the four slices in apparent magnitude, the 
amplitude of w{9) measured in the Subaru data agrees 
well with the measurements in the ACS. At very small 
angular separations the ACS data are higher than the 
ground-based results. For the two faintest bins (23 < 
i < 24 and 24 < i < 25) the agreement between the ACS 
measurements and the Subaru measurements is excellent. 

We now examine more closely the galaxy correlation 
function measured from the ACS catalogs. As we have 
already seen, thanks to the excellent resolution of the 
ACS images, we are able to measure clustering ampli- 
tudes to small separations, on the order of 1". Con- 
versely, as a consequence of the large areal coverage of 
the ACS COSMOS field, we can also measure amplitudes 
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Fig. 1. — Galaxy counts for the three catalogues presented in this 
paper: ACS (solid line); CFHTLS-T03 (open squares) and Subaru 
(red filled triangles) . For ref erence, we also show galaxy counts 
extracted from the HDF-N bv lMetcalfe et al.l <200ll) 



0.1 r 

0.01 

0.001 
1 

0.1 
0.01 k 



0.001 



"'I I 



tr i i i q 




21 < i < 22 
llllll| I I lll| I lllll 
23 < i < 24 




22 < i < 23 

| iiiii| i n iiiii| i i n iiii| i In4 

24 < i < 25 




0.001 0.01 0.1 1 0.001 0.01 0.1 1 
0(o) 0(0) 



Fig. 2. — The angular correlation function, w(8) as a func- 
tion of angular separation, 9 in four slices of apparent magni- 
tude. In each panel we show three different measurements from 
three different catalogs: ACS (connected lines); Subaru (triangles); 
CFHTLS/Megacam (stars). 

to large angular separations. In Figure[3]we show the an- 
gular correlation function w(0) as a function of angular 
separation for bright and faint samples with 21 < i < 22 
and 24 < i < 25 respectively. The dotted line shows 
the best-fitting lines with an integral constraint correc- 
tion applied. For the bright bin, we find a best fitting 
slope of -0.59 ± 0.05; for the fainter bin, -0.47 ± 0.02. 
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Fig. 3. — The amplitude of the angular correlation function w, as 
a function of angular separation, 8, for the ACS catalog. Measure- 
ments for galaxies selected in the magnitude ranges 21 < i < 22 
and 24 < i < 25 and are presented. The dotted lines show the 
best-fitting power-law correlations (with slopes —0.59 and —0.47, 
respectively) and the integral correction for the COSMOS field 
included. The dashed line is plotted at 1'. There seems little evi- 
dence here for any deviation from power-law correlations, although 
the brighter bin does hint at an inflection around 1 arcmin. The 
correlation function for the fainter bin is clearly flatter. 

We find unambiguously that in this magnitude limited 
sample, the galaxy correlation function becomes flatter 
towards fain ter magnitude bins, in agreement with pre- 
vious works (|McCracken et al"]|2001[ ). 

4. COMPARISONS WITH SIMULATIONS 

In the previous Sections we established that measure- 
ments between the different datascts arc consistent. In 
this Section, we compare our me asurements to those 
made on catalogs produced by iKitzbichler &: W hite 
(200$. These catalogs were created using a semianalytic 
model to simulate galaxy formation within the evolv- 
ing halo population of the extremely large Millennium 
Simulation. These simulated Universes were then 'ob- 
served' to produce light-cones, which can be then used 
to produce obs ervations of identica l geometry to those of 
real c atalogs. dBlaizot et alj|2005t IKitzbichler k White! 
I200H . Each COSMOS light-cone covers 2 deg 2 . Figureg] 
shows, as before, the amplitude of w{&) as a function of 
angular separation in four magnitude slices. Points with 
error bars show measurements from the ACS data. The 
solid line shows the average of measurements made from 
twenty light-cones extracted from the millennium simu- 
lation; the dotted lines show the amplitude of the ±ltr 
error bars. 

For the two brighter slices, 21 < i < 22 and 22 < i < 
23, the agreement between the simulations and obser- 
vations at intermediate and small scales is remarkably 
good. At fainter magnitudes, and at larger scales, the 
amplitudes measured in the COSMOS catalogs are con- 
sistently higher than the prediction of the simulations. 
Since the simulations allow us to assess the amplitude of 



Fig. 4. — The amplitude of the angular correlation function 
■w, as a function of angular separation, 8, for the ACS catalog. 
The lines show measurements made using mock catalogs extracted 
from the millennium simulation. The dotted lines show the one 
sigma scatter of the results for individual mock catalogs around 
their mean computed from the variance over twenty mock catalogs. 
Note that these mock results automatically incorporate integral 
corrections due to the finite field size. 

cosmic variance directly, we can easily conclude that this 
discrepancy cannot plausibly be attributed to the COS- 
MOS field having above-average clustering. It therefore 
appears that the model predictions for w(0) at i ~ 24 
are too low; we discuss this further below. 

5. COMPARISONS WITH LITERATURE MEASUREMENTS 

From our measurements it is clear that the form of 
w(8) in the COSMOS field does not correspond to a sim- 
ple power law with a slope that is independent of median 
magnitude of the sample. In the past, determinations 
of galaxy clustering amplitude were usually given at a 
fixed angular scale over a small range of angular separa- 
tions. From our data, it is clear that fitting w(6) only 
at small angular separations will result in different fitted 
amplitudes as compared to a fit over the entire range. 
Therefore, in order to compare with results presented in 
the literature, we choose to carry out fits over a similar 
range of angular separations. 

In Figure [5] we show the fitted amplitudes of w(9) as 
a function of the median magnitude of each slice. In 
comparing with literature measurements, we see that at 
bright magnitudes (i ~ 20) our measurements are ap- 
proximately in agreement with the values presented in 
the literature. However, by i ~ 23, the COSMOS field 
measurements are significantl y higher than, for exam- 
ple the CFDF measurements of lMcCracken et al.l (|2001h . 
This is true at all angles and for all measurements, 
whether they are from the ground based or space-based 
catalogs. 

Photometric r edshifts are available fo r all objects in the 
COSMOS field (jMobasher et al.ll200l . The slice with 
22 < i < 23 has a median redshift of z ~ 0.8; the next 
faintest slice is at z ~ 0.9. If we examine Figure 2 of 
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Fig. 5. — Comparisons of the COSMOS results for the depen- 
dence of the amplitude of w(6) on depth, compared with existing 
measurements and with the mock results. We choose to specify the 
amplitude at 1 arcmin, based on a <5 = —0.8 power-law fit around 
this point. The results are highly insensitive to this choice of slope. 



iScoville et all (|2007bD we can see there is a significant 
over-density in the redshift range 0.7 < z < 0.9. It is 
possible that this structure could contribute to the en- 
hanced signal on large scales seen in our data in these 
slices. 

Some independent evidence exists that this structure 
raises the COSMOS clustering amplitude above the en- 
semble average. In Figure [5] we have plotted the aver- 
age amplitude of id as a function of apparent median 
magnitude for the four independent deep fields of the 
Canada-France Legacy Survey (one of which is the D2 
field described above), totalling an effective area of 3.2 
deg 2 . Error bars correspond to the variance over all four 
fields. Interestingly, at bright magnitudes, the CFHTLS 
magnitudes agree with the COSMOS measurements, and 
other literature values; at intermediate magnitudes (22 < 
i < 23) they are between the COSMOS values and those 
from other surveys, whereas at fainter magnitudes they 
agree perfectly with the COSMOS measurements, pre- 
sumably because the median redshift probed by both 
surveys at these magnitude limits is beyond that of the 
rich structure in the COSMOS field. The size of the cos- 
mic variance error bars on the CFHTLS measurements 
are also largest at intermediate magnitudes. These mea- 
surements are discussed in greater detail in a forthcoming 
paper (McCracken et al., in preparation). 

6. DISCUSSION AND CONCLUSIONS 

This paper has presented measurements of the angu- 
lar two-point correlation function in the COSMOS field, 
w(9), and how it depends on i-band magnitude depth. 



We have shown that consistent results are obtained using 
three independent datasets: HST ACS; Subaru Suprime- 
Cam; and CFHT MegaCam. The results agree well at 
bright magnitudes (i ~ 22) with previous measurements 
and with the predictions of semi-analytic mock catalogs 
constructed from the Millennium Simulation (MS). The 
only slight caveat here is that the results at 9 > 10 ar- 
cmin are at the high end of the MS predictions, which 
may reflect a single rich z ~ 0.8 structure in the field. 

At fainter magnitudes, however, a different picture 
emerges. By the time we reach the 24 < i < 25 bin, 
the COSMOS measurements are consistently a factor 2 
higher than the MS predictions at all angles. Moreover, 
the COSMOS measurements are consistent with the four- 
field average of measurements from the CFHTLS survey. 
This discrepancy is well beyond the compass of cosmic 
variance from limited numbers of rich structures, as mea- 
sured via the ensemble of simulations. Thus, barring 
some undetected systematic that is consistent between 
all the datasets we have used, the conclusion must be 
that the MS predictions are too low at these magnitude 
levels. This could arise in a number of ways: the pre- 
dicted degree of bias at high redshifts might be too low; 
the luminosity function might be incorrect, resulting in 
too high a predicted mean redshift at these depths; or 
alternatively the MS may miss foreground pairs of intrin- 
sically faint galaxies because of the resolution limit of the 
simulation. The first possibility is particularly interest- 
ing given the current debate over the normalization of the 
primordial power spectrum, <rg. The MS used erg = 0.9, 
whereas WMAP favors a smaller result - perhaps as low 
as as — 0.7 (Spergel et al. 2006). Since high-redshift 
galaxies are strongly biased already, a reduced erg will 
in fact boost the predicted galaxy clustering (for a given 
galaxy mass) . These issues will be explored further in fu- 
ture papers, where we make direct use of the photometric 
redshift data in the COSMOS field. 
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